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We present the status of the description of hadron production in central 
nucleus-nucleus collisions within the statistical model . Extending previous 
studies by inclusion of very high-mass resonances (m> 2 GeV) , and the up- 
to-now neglected scalar a meson leads to an improved description of the 
data. In particular, the hitherto poorly reproduced energy dependence of 
the K + /ir + ratio at SPS energies ( "the horn" ) is now well described through 
the connection to the hadronic mass spectrum and, implicitly, Hagedorn's 
limiting temperature. 

PACS numbers: 25.75.Dw,25.75.Nq 

1. Introduction 

One of the major goals of ultrarelativistic nuclear collision studies is to 
obtain information on the QCD phase diagram pQ. A promising approach 
is the investigation of hadron production. Hadron yields measured in cen- 
tral heavy ion collisions from AGS up to RHIC energies can be described 
very well [21[3l[ll[5l[6l[7l[8l[9l [TOl HH [12l [131 HU within a hadro-chemical 
equilibrium model. In our approach [21 [U [12] the only parameters are 
the chemical freeze-out temperature T and the baryo-chemical potential /if, 
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(and the fireball volume V, in case yields rather than ratios of yields are 
fitted) ; for a review see [15] . 

The main result of these investigations was that the extracted temper- 
ature values rise rather sharply from low energies on towards ^/snn —10 
GeV and reach afterwards constant values near T=160 MeV, while the bary- 
ochemical potential decreases smoothly as a function of energy. The limiting 
temperature [Hj behavior suggests a connection to the phase boundary and 
it was, indeed, argued [T7| that the quark-hadron phase transition drives the 
equilibration dynamically, at least for SPS energies and above. Considering 
also the results obtained for elementary collisions, where similar analyses 
of hadron multiplicities, albeit with several additional, non-statistical pa- 
rameters (see [181 H] an d refs. therein), yield also temperature values in 
the range of 160 MeV, alternative interpretations were put forward. These 
include conjectures that the thermodynamical state is not reached by dy- 
namical equilibration among constituents but rather is a generic fingerprint 
of hadronization |20| [21J, or is a feature of the excited QCD vacuum |22j . 
The results presented below lend further support to the interpretation that 
the phase boundary is reflected in features of the hadron yields. 

While in general all hadron yields are described rather quantitatively, a 
notable exception was up-to-now the energy dependence of the K + /ir + ratio 
which exhibits a rather marked maximum, "the horn" [25], near y / ijv r /v ~ 10 
GeV [23J. Predicted first within a model of quark-gluon plasma (QGP) for- 
mation [25] , the existence of such a maximum was also predicted [23] within 
the framework of the statistical model, but the observed rather sharp struc- 
ture could not be reproduced [12]. Other attempts to describe the energy 
dependence of the K + /ir + ratio within the thermal model [111114] were also 
not successful, except when an energy-dependent light quark fugacity was 
used as an additional parameter [14]. Furthermore, all attempts to repro- 
duce this structure within the framework of hadronic cascade models also 
failed, as is discussed in detail in [23] . As a consequence, the horn structure 
is taken in [23] as experimental evidence for the onset of deconfinement and 
QGP formation, and as support for the predictions of [25 j . We have recently 
shown [30] that, employing an improved hadronic mass spectrum, in which 
the a meson ]29] and many higher-lying resonances are included, leads to 
an increase of about 16% for the calculated pion yields. This increase levels 
off near the point where the temperature reaches its limiting value, thereby 
sharpening the structure in the K + /ir + ratio, as will be shown below. For 
the A hyperons, the new high mass resonances lead to an increase in the 
calculated production of about 22%. An increase of up to 6% is observed 
for protons, while for kaons this increase is about 7%. 

In Fig. [1] we present a comparison of measured and calculated hadron 
yields at the energies of ^/saw=7.6 GeV (beam energy of 30 AGeV at SPS) 



therm printed on April 8, 2009 



3 



10 









\&=7-6GeV _ z 


















































[ 






1 


r 




, ■ Data (NA49) 






- Model, ]( Z /N =16.3/8 




■ 


: t= 


138 MeV.n =380 MeV, V= 


1380fm 3 



w 



10 



Data 



STAR 

PHENIX 

BRAHMS 

Mbde),x s /N,,|=S 
T=164 MeV.n 



s MNI =200 GeV 



tr 



30MeV, V=1950 fm 3 



r 



%* % K + K p 



A A E 



IfflCKppAAHSfltcId K* E* A* 



Fig. 1. Experimental hadron yields and model calculations for the parameters of 
the best fit at the energies of 7.6 (left panel) and 200 GeV (right panel; the £1 yield 
includes both Vl~ and f2 + ). 

and y / ijv r /v=200 GeV. The model is successful in reproducing the measure- 
ments and this applies to all energies, from 2 AGeV beam energy (fixed 
target) up to the top RHIC energy of y / sjv r /v=200 GeV. The reduced \ 2 
values are reasonable. In most cases the fit quality is improved compared 
to our earlier analysis [T2], even though the experimental errors are now 
smaller. A disagreement between the experiments is seen at the top RHIC 
energy for pions and protons, see Fig. (H which is the reason of the large 
reduced x 2 - A fit of ratios is in this case more suited, but we note that a 
fit of the STAR yields alone gives T=162 MeV, fi b =32 MeV, V=2400 fm 3 , 
with a very good x 2 /Ad/=9.0/ll. The resonances were not included in the 
fits, but are quite well reproduced by the model. 

An important result of our analysis is that the resulting thermal param- 
eters are close to those obtained earlier |12j and are in agreement with other 
recent studies |13|. [28] . This indeed confirms that the common practice of 
including in the thermal codes hadrons up to masses of 2 GeV (for instance 
in the publicly-available code THERMUS [31] ) does not lead to significantly 
biased fit parameters. Nevertheless, there are small variations. In Fig. [2] 
we present the energy dependence of T and ^ in comparison to our ear- 
lier results [T2]. We have parametrized T as a function of y / J/v r /v with the 
following expression!]]: 



T = T U 



1 1 + exp(2.60 - ln(VsAw(GeV))/0.45) 



(1) 



For Hb, there is no need to change our earlier [T2] parametrization: /iiJMeV] = 

1 inr> 



1 + 0.286-y/sjvjvtGoV) 
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with the "limiting" temperature T/j m =164 MeV. This value is slightly higher 
compared to our earlier value of 161±4 MeV |12j due to the higher tem- 
peratures presently derived for the RHIC energies. The approach to T^ m is 
presently more gradual compared to our earlier parametrization. 




Vs NN (GeV) Vs NN (GeV) 

Fig. 2. The energy dependence of the temperature and baryon chemical potential 
at chemical freeze-out. The results of the new fits ; 30| are compared to the values 
obtained in our earlier study [12 . The lines are parametrizations for T and Hb (see 
text). 

The values of /ib extracted for the two lowest SPS energies deviate some- 
what from the continuous trend suggested by all the other points. At these 
energies the fit to data in full phase space does lead, as expected, to larger 
values of m,, which do fit in the systematics. At 40 AGeV (y / J/vw r =8.8 GeV) 
the resulting values of T and from the fit of 4ir data are very similar to 
those obtained from midrapidity data. 

The volume at chemical freeze-out (corresponding to a slice of one unit 
of rapidity, dV/dy) is shown in Fig. as a function of energy. The values 
extracted directly from the fits of particle yields are compared to the values 
obtained by dividing the measured charged particle rapidity densities with 
densities calculated employing the above parametrizations of T and fib- As 
expected, the two methods give very similar results. The chemical freeze- 
out volume is compared to the kinetic freeze-out volume extracted from 
Hanbury Brown and Twiss (HBT) measurements, Vhbt |32j . Note that, to 
relate quantitatively the magnitude of the rapidity density of the chemical 
freeze-out volume dV/dy to the volume at kinetic freeze-out determined 



therm printed on April 8, 2009 



■5 



E 

aT4000 
E 

£ 3500 
3000 
2500 
2000 
1500 
1000 
500 



10 10 2 

Vs NN (GeV) 

Fig. 3. The fireball volume at chemical freeze-out compared to the volume extracted 
from the HBT measurements [3"2] . 

from HBT measurements, one needs to map rapidity onto space (see the 
discussion in [32J ) . Here we only observe that the energy dependence of the 
two observables exhibit a similar non-monotonic behavior. It is interesting 
to note that the minimum of dV/ dy occurs at the ^/snn corresponding to 
the saturation of T. 

We employ the above parametrizations of T and /if, to investigate the 
energy dependence of the production yields of various hadrons relative to 
pions, shown in Fig. HJ In particular, the K + /tt + ratio shows a rather pro- 
nounced maximum at a beam energy of 30 AGeV [23J, and the data are well 
reproduced by the model calculations. In the thermal model this maximum 
occurs naturally at ^snn —8 GeV [23]. It is due to the counteracting effects 
of the steep rise and saturation of T and the strong monotonous decrease 
in //ft. The competing effects are most prominently reflected in the energy 
dependence of the A hyperon to pion ratio (right panel of Fig. (3|), which 
shows a pronounced maximum at y / Siv r /v —5 GeV. This is reflected in the 
K + /ir + ratio somewhat less directly; it appears mainly as a consequence of 
strangeness neutrality, assumed in our calculations. 

The model describes the K + /tt + data very well over the full energy 
range, as a consequence of the inclusion in the code of the high-mass reso- 
nances and of the a meson, while our earlier calculations [12] were overpre- 
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Fig. 4. Energy dependence of hadron yields relative to pions. 



dieting the SPS data. At RHIC energies, the quality of the present fits is 
essentially unchanged compared to [12] , as also the data have changed some- 
what. The model also describes accurately the A/tt~ measurements as well 
as those for other hyperons. We note that the maxima in the various pro- 
duction ratios are located at different energies [33]. The model calculations 
reproduce this feature in detail. 

The calculated K + /ir + ratio is likely to decrease further at energies 
beyond the maximum and the peak is likely to sharpen somewhat [30J if our 
presumably still incomplete knowledge of the hadronic spectrum for masses 
larger than 2 GeV would improve. The uncertainty of the calculations due to 
the mass and width of the a meson are at the percent level only. Another 
few percent uncertainty, which is difficult to assess quantitatively, arises 
from the unknown branching ratios of the high-mass resonances. 
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In summary, our results [30] demonstrate that by inclusion of the a me- 
son and many higher mass resonances into the resonance spectrum employed 
in the statistical model calculations an improved description is obtained of 
hadron production in central nucleus-nucleus collisions at ultra-relativistic 
energies. A dramatic improvement is visible for the K + /tt + ratio, which 
is now well described at all energies. The "horn" finds herewith a natural 
explanation which is, however, deeply rooted in and connected to detailed 
features of the hadronic mass spectrum which leads to a limiting temper- 
ature and contains the QCD phase transition [16]. Our results strongly 
imply that hadronic observables near and above the horn structure at a 
beam energy of 30 AGeV provide a link to the QCD phase transition. Open 
questions are whether the chemical freeze-out curve below the horn energy 
actually traces the QCD phase boundary at large values of chemical po- 
tential or whether chemical freeze-out in this energy range is influenced by 
exotic new phases such as have been predicted in [5S] . In any case these are 
exciting prospects for new physics to be explored at RHIC low energy runs 
and, in particular, at the high luminosity FAIR facility. At the high en- 
ergy frontier, the measurements at LHC will be a crucial test of the present 
picture. 
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